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FOREWORD 

This  report  was  prepared  by  A.  6 Gerardl,  Aerospace  Engineer  in 
the  Loads  and  Response  Prediction  Group  of  the  Structural  Mechanics 
Division  of  the  Air  Force  Flight  Dynamics  Laboratory  at  Wright-Patterson 
Air  Force  Base,  Ohio.  The  work  described  herein  Is  a part  of  the  Air 
Force  Systems  Command  exploratory  development  program  to  predict 
aircraft  dynamic  loads  during  ground  operations.  The  work  was  directed 
under  Project  1367,  "Structural  Integrity  for  Military  Aerospace 
Vehicles,"  Task  136701,  "Structural  Flight  Loads  Data." 

This  report  covers  work  done  In  the  period  from  September  1975  to 
August  1976. 


*!:CESStON  f(f 

KTIS  Vette 

OOC  fkiblMOM  o 

UfMMKOUKCtO 
MTIflCATkMI 


i?r 


nwmwuindB 


A 


ill 


AFFOL-TR- 77-37 


TABLE  OF  CONTENTS 


Section 


I INTRODUCTION 

1.  Purpose  of  the  Study 

II  MATHEMATICAL  MODEL 

1.  General  Airplane/Runway  Model 

2.  Rigid  Body  Equations  of  Motion 

3.  Flexibility  Equations  of  Motion 
4=  Solution  Technique 

III  COMPUTER  PROGRAM 

1 . Output  Format 

IV  DISCUSSION  OF  SIMULATKWS 


SUMMARY  AND  CONCLUSIONS 


APPENDIX 

A DEVELOPMENT  OF  EQUATIONS  OF  MOTION 
B LISTING  OF  COMPUTER  PROGRAM  TAX2 

C LISTING  OF  AIRPLANE  DATA 

D LISTING  OF  RUNWAY  PROFILE  DATA 
BIBLIOGRAPHY 
REFERENCES 


PASS  PiilM 


iiiitiaiiiii*iiitMiiiMi!iiaiMiiMiiM<iiiiiiiiiaiiiaB4wiM^^^^  * 


AFFDL-TR-77-37 


LIST  OF  ILLUSTRATIONS 

Figure 

Page 

1. 

Accepted  Military  Human  Tolerance  Vertical 

Vibration  Criterion 

2 

2. 

Typical  Single  Acting  Oleo  Pneutnatic  Landing  Gear 

Strut 

4 

3. 

Source  Deck  Setup 

16 

4. 

Typical  Cal comp  Plotted  Output 

18 

5. 

Boeing  727-100  Constant  Speed  Taxi  Simulation  over 
the  Washington  National  Profile  Without  a Roll 

Degree  of  Freedom 

22 

6. 

Boeing  727*100  Constant  Speed  Taxi  Simulation  over 
the  Washington  National  Profile  With  a Roll  Degree 
of  Freedom 

23 

7. 

PSD  of  Washington  National  Airport  Runway  36 

24 

It-, 

Boeing  727-100  Traversing  a (1-cos)  dip  head-on 

25 

9, 

bv-eind' 727-100  Traversing  a (l-cos)  dip  at  a 

45*  angle 

26 

10. 

C-9A  Traversing  a (1-cos)  dip  at  a 45*  angle 

28 

11. 

AMST  Traversing  a (1-cos)  dip  at  a 45*  angle 

29 

12. 

Boeing  727-100  Taking  Off  from  Washington 

National  Airport  With  the  Roll  Degree  of 

Freedom  Included 

30 

13. 

AMST  Taking  Off  from  Washington  National 

Airport  With  the  Roll  Degree  of  Freedom  Included 

31 

14. 

C-9A  Taking  Off  from  Washington  National  Airport 

With  the  Roll  Degree  of  Freedom  Included 

32 

15. 

PSD  of  Washington  National  Runway  36  and  two 

Typically  Smooth  Runways 

33 

16. 

Measured  Response  of  a Boeing  727-100  Takeoff  at 
Washington  National  Airport  Runway  36 

34 

17. 

C-9A  with  Flexible  Wings  Taxiing  over  a (1-cos) 
dip  at  a 45°  angle 

36 

v1 


AFFDL-TR.77-37 


LIST  OF  ILLUSTRATIONS  l! Conti r4jed) 

Figure  page 

18.  C-9A  with  Flexible  Wings  Taking  Off  from 

Washington  National  Runway  36  37 

A-1.  Description  of  Asynwetrical  Mathematical  Model  42 

C-1.  Three  View  of  Boeing  727-100  63 

C-2.  Three  View  of  McDonnell -Douglas  C-9A  $6 

0-1.  Washington  National  Runway  36  with  Linear 

Trend  Removed  72 


vii 


*,IST  or  TABLES 


AfFOL*TB-77-3? 


Table 

1 Descrlptioii  of  Input  Aircraft  Data 

2 Typical  CoetM/<»»r  Output  Listing 
St«aiy  of  Sl^latlons 

Co^rfsons  of  Slmulate.<l  and  ExperlEcntal  Data 


1 


AFFDL-TR-77-37 


SECTION  I 
INTRODUCTION 

A coBwon  problem  that  can  occur  during  takeoff  and  taxiing  operations 
of  aircraft  is  high  acceleration  levels  caused  by  a rough  runway.  Due 
to  these  accelerations,  runway,  must  be  evaluated  with  respect  to  roughness 
in  order  to  ensure  timely  pavemer.t  maintenance  to  control  aircraft 
structural  loads  and  f'a-.igue.  Also,  rough  runways  adversely  affect  the 
ability  of  the  crew  members  by  reducing  instrument  readability  and 
crew  comfort.  Figure  1 show:  the  current  criterion  (Reference  1)  used 
to  set  maximun  allowable  vertical  acceleration  levels  from  a human 
comfort  standpoint.  Reference  2 addressed  the  runway  roughness  problem 
at  considerable  length  and  contains  the  devclo^ent  of  a mathematical 
model  and  subsequent  computer  program  called  ”VAXP  to  simulate  the 
dynamic  response  of  military  aircraft  to  runway  roughness  on  a synmetrical 
runway.  For  a symmetric  runway,  only  one  runway  profile  is  required. 
Normally  this  is  sufficient  for  representing  a paved  runway.  With 
the  advent  of  the  AMST  (Advanced  Medium  STOL  Transport)  and  is  some  cases 
with  conventional  airplanes  operating  off  of  semiprepared  or  very  rough 
paved  surfaces,  the  rolling  motion  of  an  aircraft  became  significant. 

This  rolling  motion  was  the  result  of  operating  the  aircraft  on  an 
asymmetric  runway.  Therefore,  in  order  to  properly  simulate  this 
response  it  became  necessary  to  include  the  runway  profile  encountered 
by  each  landing  gear. 


1 . PURPOSE  OF  THIS  STUDY 


The  purpose  of  this  study  is  to  develop  a computer  program,  capable  i | 
of  simulating  an  aircraft  during  constant  speed  taxi  or  takeoff  from  | 
ninways  that  are  asywoetrical . j 
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SECTION  II 
MATHEMATICAL  MODEL 

The  airplane/ runway  mathematical  model  used  for  this  study  was  the 
basic  mathematical  model  developed  in  Reference  2.  A detailed  description 
of  the  components  that  make  up  this  general  model,  as  well  as  the 
assumptions  made  are  shown  in  Reference  2.  This  report  presents,  in 
simmary  form,  the  landing  gear  strut  and  tire  representation,  the 
airplane  rigid  body  and  flexible  body  representation,  the  runway 
profile  representation,  the  equations  of  motion,  and  the  solution 
technique. 

1.  GENERAL  AIRPLANE/RUNWAY  MODEL 

The  general  model  represents  an  asymnetrical  body  with  a nose  gear 
and  a right  and,  left  main  landing  gear.  Each  landing  gear  strut  is 
assumed  to  have  point  contact  with  the  profile  and  it  is  assumed  that 
each  landing  gear  traverses  a different  profile.  Aerodynamic  lift  and 
drag  are  modeled,  and  thrust  is  applied  at  the  aircraft's  center  of 
gravi ty. 

The  airplane  is  free  to  roll,  pitch,  plunge,  and  translate  horizontally 
down  the  runway  and  each  landing  gear  unsprung  mass  is  free  to  translate 
vertically.  To  these  rigid  body  degrees  of  freedom,  up  to  30  flexible 
modes  of  vibration  are  included.  This  airplane  motion  is  controlled 
by  the  landing  gear  strut  forces,  lift,  drag,  thrust,  and  the  resisting 
parameters  of  aircraft  mass  and  inertia. 

The  landing  gear  struts  are  nonlinear,  single  acting  oleo  pneumatic 
energy  absorbing  devices  {Figure  2)  and  are  represented  in  the  model 
as  the  sum  of  the  three  forces;  pneianatic,  hydraulic,  and  strut  bearing 
friction  forces.  The  pneumatic  force,  which  is  the  largest  of  the  three 
is  represented  by  the  equation: 
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Figure  2.  Typical  Single  Acting  01 eo  Pneumatic 
Landing  Gear  Strut 
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where: 

P = fully  extended  strut  pressure 
V = fully  extended  strut  volume 
A = pneumatic  piston  area 
S = strut  stroke 

The  damping  force  is  given  by  the  equation: 

^ 2 (C^  A^)^- 


(2) 


where: 

= density  of  the  hydraulic  fluid 
Ajj  = the  hydraulic  piston  area 

A = effective  orifice  area  (constant  orifice  minus  metering 
“ pin  area) 

Cjj  * orifice  coefficient  (use  0.9) 

S = strut  piston  velocity 

The  third  strut  force  is  the  strut  bearing  friction  force  and  is 
neglected  in  the  model  because  the  force  is  small  for  symmetrically 
loaded  struts.  (See  Reference  2). 


The  tire  force  is  represented  by  the  linear  equation: 

» k T 
T D 

where; 

Tp  = tire  deflection 
k * linear  tire  spring  constant 


The  runway  elevation  data  is  input  into  the  model  in  two  foot 
Increments.  The  profile  is  made  continuous  by  fitting  the  following 
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polynomlnal  through  the  three  elevation  data  points  a.id  the  slope  at 
the  end  of  the  previous  profile  segment: 

y{x)  - + a2X  + a,x^  + a^x^  (4) 

where*. 

9 -i  A ~ coefficients  derived  from  the  elevation  and 
**^'^*‘*  slope  data 

This  is  done  for  each  of  the  three  lines  of  runway  profile  data. 

2.  RIGID  BODY  EQUATIONS  OF  MOTION 

The  differential  equations  of  motion  for  the  mathematical  model 
were  derived  by  application  of  the  Lagrange  equations  (See  Appendix  A). 
The  general  form  of  these  equations  is  shown  below  and  corresponds  to 
the  notation  shown  in  Figure  A-1  in  Appendix  A. 

Z » (Fgj^  + Fg2  *'^83  + L - Ic.g.»  vertical  acceleration]  (5) 

*«  (F^j^  - Fgj  - (imsprung  mass  vertical  acceleration]  (6) 

^2  “ ^'^tZ  “ ^s2  ” ^2^^^  (unsprung  mass  vertical  acceleration]  (7) 

•• 

Zj  » (F^„  - Fgj  - Wjl/Mj  (unsprung  mass  vertical  acceleration]  (8) 

0 * ^^sl^  ^s2®  ^Td'^I  ” 'sS^^^^yy  (pitching  acceleration]  (9) 

** 

* (F^2  - (rolling  acceleration]  (10) 

X ■ (Fj  - F.JP  - F^)/(M^g)  (horizontal  translation  acceleration]  (11) 
where: 

^sl*  ^s2»  ’^sS  > total  landing  gear  strut  forces 
^ti*  ^c2»  ^t3  “ forces 

M 1 W,  I . I « aircraft  mass,  weight*  and  pitching  and 
yy  XX  inertias 

*The  subscript  1,  2 and  3 corresponds  to  the  nose,  right  main  and 
left  main  landing  gears  respectively. 
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u u u 
1’  2'  3 

A.  B,  C,  £, 


upsprung  landing  gear  weights 
noisent  arms 


L,  F__,  F._  « lift,  thrust,  and  tire  and  aerodynamic 

1 * 1/  • aL*  • 

drag  forces 

(F^  and  act  through  the  center  of  gravity] 


3.  FLEXIBILITY  EQUATIONS  OF  MOTION 
"i'^i  “ ^il^sl  * ^i2^s-2  ^13^83  ’ ^^"l^i^i  ■ “i^^i**! 


where: 


M.  = the  generalized  mass 

^11  ‘ ^12*  ^13  " deflections  at  gear  location  1,  2 and  3 


0)^  * modal  frequency 


5 = damping  factor 

Ql  = Ri  = generalized  coordinates  and  their  time  derivatives. 
The  sign  convention  Is  as  follows: 


z * Vertical  Displacement 

4 up 

0 “ Pitch 

+ nose  down 

^ = Roll 

4-  roll  right 

q » Deflection  Dae  to  Bending 

4-  up 

X Horizontal  Translation 

4-  forward 

4.  SOLUTION  TECHNIQUE 

The  technique  used  for  solving  the  coupled  nonlinear  differential 
equations  of  motion  that  describe  the  simulated  aircraft  is  a three- 
term  Taylor  series.  For  example,  the  equation: 

X ■ -cx  -kx 

The  three  term  Taylor  series  representations  can  be  written  as 

*(M)  ■ ’‘(I)  * *<I)  <“>  ♦ *(I)  ^ <“> 

where:  1 = 1-^ 
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The  values  for  x,  x and  x from  the  previous  step  are  substituted  into 
Equation  13  and  a new  value  for  x is  obtained.  Differentiating 
Equation  13  we  obtain  for  the  velocity  x,  the  expression*. 


’'(i+l) 


’‘(I)  *(!) 


The  values  for  x and  x are  then  substituted  into  Equation  14  and 
a new  value  of  x is  found.  This  entire  process  is  repeated  with  the 
new  values  of  x and  x to  obtain  the  next  point  in  the  solution. 
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SECTION  m 
COMPUTER  PROGRAM 

The  computer  program,  TAX2,  vrtiich  computes  the  dynamic  response 
of  a flexible  aircraft  to  an  asymmetrical  runway  profile,  consists  of 
one  basic  program  and  several  subroutines.  A complete  listing  of  the 
program  is  contained  in  Appendix  B.  Table  1 contains  a description  of 
the  aircraft  input  data  and  Figure  3 shows  the  source  deck  setup  for 
use  on  the  CDC  6600  computer  at  Wright-Patterson  AFB,  Ohio. 


1 . OUTPUT  FORMAT 

The  results  of  the  calculations  are  presented  as  both  a printed 
output  and  1 time  history  plot.  The  printed  output  lists  the  value 
of  fifteen  parameters  each  0.01  second.  A sample  of  this  listed  output 
is  shown  in  Table  2.  The  fifteen  parameters  listed  in  the  heading  are; 

XMAINL  - Left  Main  landing  gear  strut  deflection  (inches) 

XMAINR  - Right  Main  landing  gear  strut  deflection  (inches) 

XNOSE  - Nose  gear  strut  deflection  (inches) 

ZPW.  - Left  Main  landing  gear  runway  elevation  (inches) 

ZPMR  - Right  Main  landing  gear  runway  elevation  (inches) 

ZPN  - Nose  landing  gear  runway  elevation  (inches) 

BETADD  - Rolling  acceleration  (^)  (rad/sec^) 

THETAOD  - Pitching  acceleration  (o)  (rad/sec^) 

BETA  - Roll  angle  (4»)  (rad) 

THETA  - Pitch  angle  (o)  (rad) 

SPEED  - Aircraft  velocity  (ft/sec) 

DISTANCE  - Distance  traveled  down  the  runway  (feet) 

TIME  - Real  time  (seconds) 

C6ACC  - Center  of  Gravity  Vertical  Acceleration  (g‘s) 

PSA  - Pilot's  Station  Vertical  Acceleration  (g's) 


Figure  4 shows  a photographic  reduction  of  a typical  Calcomp- 
plotted  time  history.  ”ni1s  figure  depicts  a Boeing  727-100  tar ling 
at  50  fps  over  a l-cos  bump  at  a 45^  anglt^  of  approach.  The  plotted 
output  includes  titles  showing  the  airplane  simulated,  its  gross  weight. 
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TABLE  1 

DESCRIPTION  OF  INPUT  AIRCRAFT  DATA 
Section  1 (cards  1-5)  - General  Airplane  Data 


Card 

Column 

Format 

Variable 

Name 

Data  for 
McDonnell- 
Oouglas  C-9A 

Definition 

Card  1 

1-80 

8A10 

PLANE 

McOonnell- 
Douglas  C-9A 

Airplane  Being  Simulated 
and  Gross  Weight 

Card  2 

1-10 

FlO.l 

W 

108000. 

Vehicle  Weight  (lbs) 

11-20 

FlO.l 

A 

51.6 

Distance  Main  Gear  to  CG 
(in) 

21-30 

FlO.l 

B 

589.4 

Distance  Nose  Gear  to  CG 
(in) 

31-40 

FlO.l 

MMI 

20800000. 

Pitch  Momeot  of  Inertia 
(lb  in  sec^) 

41-52 

FI  2.0 

WS 

96. 

Wing  Station  of  Main  Gear 
(in) 

53-64 

F12.0 

MHIR 

8000000. 

Roll  Moment  of  Inertia 
(lb  in  sec^) 

Card  3 

1-10 

11-20 

F10.2 

F10.2 

PSARM 

TAILRM 

607.0 

318.5 

Distance  of  Pilot  Station 
to  CG  (In) 

Distance  of  Tail  Station 
to  CG  (in) 

Card  4 

1-10 

F10.2 

SPEED 

50. 

Initial  Velocity  of  Air- 
plane (ft/sec) 

11-20 

F10.2 

THRUST 

29000. 

Total  Airplane  Thrust 
(lbs) 

21-30 

F10.2 

TAKEOFF 

285.5 

Airplane  Rotation  Speed 
(ft/sec) 
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TABLE  1 (Continued) 


Card 

Column 

Format 

Variable 

Name 

Data  for 
McOonnel 1 - 
Douglas  C-9A 

Definition 

Card  5 

1-10 

F10.4 

CL 

1.1 

Lift  Coefficient 

11-20 

FI  0.4 

AREA 

1000.7 

Wing  Area  (ft^) 

21-30 

FI  0.4 

CD 

.1 

Drag  Coefficient 

Sec 

tion  2 (cards  6-11)  - Main  and  Nose  Gear 

Card  6 

1-10 

F10.2 

WH 

957.16 

Unsprung  Weight  of  Each 

Main  Gear  (lbs) 

11-20 

F10.2 

WN 

153.43 

Unsprung  Weight  of  Nose 

Gear  (lbs) 

21-30 

F10.2 

SXM 

2. 

Number  of  Main  Gear  Struts 

31-40 

F10.2 

SXN 

1. 

Number  of  Nose  Gear  Struts 

Card  7 

1-10 

n-20 

21-30 

31-40 

F10.5 

FI  0.5 

F10.5 

F10.5 

1 1 1 I 

6.745 

8.2958 

16.5 

21.648 

Hy^aullc  Piston  Area  Nose 

Pneumatic  Piston  Area  Nose 
{1n2) 

Hydraulic  Piston  Area  Main 
(1n2) 

Pneumatic  Piston  Area  Main 

Card  8 

1-10 

F10.5 

PAOH 

120. 

Nose  Strut  Preload  Pressure 
(lbs/1n2) 

11-20 

Flo.  5 

PAON 

220. 

Mein  Strut  Preload  Pressure 

(1bs/1n2) 
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TABLE  1 (Continued) 


Data  for 


Card 

Col  u!i?n 

Format 

Variable 

Name 

McDonriell- 
Oouglas  C-9A 

Definition 

21-30 

F10.5 

VON 

126.2 

Fully  Extended  Nose  Strut 
Air  Volume  (1n3) 

31-40 

F10.5 

VCM 

366.0 

Fully  Extended  Main  Strut 
Air  Volume  (in^) 

41-50 

Flo.  5 

OAM 

.543 

Orifice  Area  Main  (In^) 

51-60 

Flo.  5 

OAN 

.442 

Orifice  Area  Nose  (1n^) 

Card  9 

1-10 

F10.3 

SLH 

85.5 

Distance  from  Axle  to  C6 
Meter) Ine 

Main  Gear  Strat  Unloaded 
(In) 

11-20 

CardJO 

FI  0.3 

SLN 

87.3 

Distance  from  Axle  to  C6 
Waterline 

Nose  Gear  Strut  Unloaded 

(1r») 

1-10 

F10.1 

7SH 

23428.6 

Main  Tire  Spring  Constant 
Per  Strut  (lbs/1n) 

11-20 

FlO.l 

TSN 

8632.5 

Nose  Tire  Spring  Constant 
Per  Strut  (Ibs/ln) 

Card  11 

1-10 

FI  0.5 

DX 

.001 

Integration  Step  Size 

Card  12 

1-5 

15 

IFPLOT 

0 

O-Plot;  1-No  Plot 

6-10 

15 

IFLIST 

0 

0-L1st;  1-No  List 

Section 

3 (cards  13-16)  - Meterlno 

Pin  Oescrlotlon 

Card  13 

1-5 

15 

NSCN 

5 

INiii^r  of  Slope  Changes 

Nose  Gear 

12 
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Card 

Column 

Format 

Variable  1 
Name  i 

* Card  14A.  14B.. 

1-10 

F10.3 

STROKN  (1) 

11-20 

F10.3 

PINDN  (1) 

Card  15 

1-5 

15 

NSCM 

* Card  16A.  16B.. 

1-10 

FI  0.3 

STROkM  (1) 

11-20 

F10.3 

PINON  (1) 

Section  4 (cards 

Card  17 

1-5 

15 

NFM 

6-10 

15 

NAFM 

*♦  Card  18A.  188.. 

♦ * * 

1-10 

F10.3 

SINAIN  (I) 

11-20 

F10.3 

3IN0SE  (I) 

21-30  no.  3 SICG  (I) 

31-40  no.3  ilTAIL  (I) 


1 (Continued) 

Data  for 
WcDonnell- 

Douglas  C-9A  Definition 


Stroke  Corresponding  to 
Metering  Pin  Diameter, 
Nose  Gear 

Metering  Pin  Diameter, 
Nose  Gear  (In) 


Number  of  Slope  Changes 
Main  Gear 


* Stroke  Corresponding  to 
Metering  Pin  Diameter, 
Nose  Gear 

* Metering  Pin  Diameter, 
Main  Gear  (In) 

17-19)  - Flexibility  Data 


7 Number  of  Sywnetrical 

Flexible  Modes 

7 Number  of  Asywaetrlcal 

Nodes 


Node  Shape  Deflection  Main 
Gear 

Mode  Shape  Deflection  Nose 

Gear 


13 


Mode  Shape  Deflection  CG 

Mode  Shape  Deflectlw 
Tall  Station 
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TABLE  1 (Continued) 
Data  for 


Card 

Column 

Variable 
Format  Name 

McDonnel 1 - 
Douglas  C-9A 

Definition 

41-50 

F10.3  SIPS 

** 

Mode  Shape  Deflection 

Pilot  Station 

**  Card 

19A.  19B 

1-15 

F15.2  6M  (I) 

** 

2 

Generalized  Mass  (lbs  sec  / 
in) 

16-25 

F10,3  OMEGA  (I) 

** 

Modal  Frequency  (rad/sec) 

**  Card  20A.  20B 

1-10 

F10.3  SILEFT  (I) 

** 

Deflection  for  Left  Main 
Gear 

11-20 

FI 0.3  SIRIGHT  (I 

j ★* ** 

Deflection  for  Right  Main 
Gear 

**  Card  21A.  21 AB..... 

1-15 

F:i5.3  GMA  (I) 

** 

Asymmetrical  Generalized 
Mass  (lbs  secvln) 

16-25 

FI 0.4  OiiEGAA  (I) 

♦* 

Asymnetrical  Modal 

Frequency  (rad/sec) 

* One  card  is  required  for  each  stroke-metering  pin  con^i nation  read 
into  the  program. 

**  One  card  is  required  for  each  flexible  mode. 

Note:  A summary  of  all  tt^ie  data  used  In  this  study  is  show:  in 
Appendix  C. 


14 


AFFOL-TR-77-37 


TABl.t  1 (Concluded) 

Runway  Profile  Magnetic  Tiipe 

The  runway  profile  is  read  into  the  program  from  a magnetic  tape 
or  permanent  file.  The  format  for  this  is  shown  below: 


Column 

Format 

Variable 

Nmne 

Definition 

Reed  1 

1-80 

8A10 

SITE** 

Runway  Profile  and  Direction 

Read  2 

1-6 

16 

NPTSS** 

Number  of  Runway  Elevation  Points 

* Read  3 

f 4« • a . . 

H*Z 

1-70 

10F7.3 

ELEVL** 

Runway  Profile  Data 

ELEV 

ELEVR 

* One  read  required  for  every  ten  runway  profile  elevation  points. 

The  process  Is  repeated  for  each  of  the  three  profiles. 

Note:  All  of  the  runway  profile  data  used  In  this  study  Is  listed 
In  Appendix  0. 
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TABLE  2 

TvPiCAL  COMPUTER  OUTPUT  LISTING 
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Plotted  Output 
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the  runway  number  and  its  location.  The  abscissa  is  the  time  axis 
annotated  every  second.  At  every  time  annotation  the  current  value  of 
aircraft  speed,  in  feet  per  second,  and  the  current  aircraft  position 
on  the  runway,  in  feet,  are  printed  out.  Runway  markers  (1,000-foot 
markers)  are  also  plotted  on  the  time  scale  to  aid  in  aircraft  positioning. 
The  plot  titled  "Nose  Gear  Track"  is  a time  history  of  the  runway 
profile  as  it  is  encountered  by  the  aircraft's  nose  gear.  The  actual 
runway  profile  is  preceded  by  100  feet  of  smooth  surface  to  allow  any 
starting  transients  to  damp  out  prior  to  encountering  the  actual  profile. 
There  are  two  aircraft  acceleration  time  histories  that  e of  particular 
interest.  One  is  the  vertical  acceleration  at  the  pilot's  station  (PSA), 
the  other  is  the  vertical  acceleration  at  the  aircraft's  center  of 
gravity  (CGA).  Each  time  history  is  banded  by  the  human  tolerance 
criterion  of  +0.4  g.  In  order  to  minimize  the  amount  of  computer 
central  memory  required  to  store  the  acceleration  time  histories, 
the  higher  frequency  components  were  effectively  filtered  out  by 
limiting  the  sampling  interval.  All  of  the  acceleration  peaks,  however, 
are  shown  on  the  plot.  It  should  be  noted  that  the  pilot  station 
acceleration  time  history  is  not  always  within  the  band  of  accepted 
human  tolerance  criteria.  Thus,  the  plot  is  very  useful  in  that  it 
provides  a graphical  record  of  the  level  of  acceleration,  and  it  shows 
which  bumps  in  the  runway  profile  caused  the  high  acceleration  response. 
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SECTION  IV 

DISCUSSION  OF  SIMUUTIONS 

Table  3 contains  a sunnary  of  the  simulations  made  In  this  study. 
Three  different  airplanes  were  simulated:  the  Boeing  727-100,  the 
McDonnell  Douglas  C-9A,  and  an  AKST  configuration.  Each  airplane  was 
simulated  traversing  two  profiles:  Washington  National  Runway  36,  and 
a l-cos  dip.  Simulations  were  made  using  mathematical  models  with  and 
without  a roll  degree  of  freedom,  I.e.  one  or  three  profiles,  and  with 
and  without  flexible  wings  so  that  comparisons  of  the  responses  could  be 
made. 

Figures  5 and  6 show  the  plotted  results  of  a Boeing  727-100 
traversing  the  Washington  National  runway  profile  without  a roll  degree 
of  freedom  and  with  a roll  degree  of  freedom  respectively.  Both  runs  were 
made  at  a constant  speed  of  125  feet  per  second,  because  this  speed  pro- 
duced higher  levels  of  vertical  acceleration  for  this  airplane.  Comparison 
of  these  two  figures  shows  a significant  Increase  In  the  vertical 
acceleration  at  the  pilot's  station  (P.S.)  while  the  aircraft  is  at 
different  locations  on  the  runway.  For  example,  at  T=46  sec.  P.S. 
acceleration  levels  more  than  doubled  when  three  lines  of  profile  were 
used.  This  Is  attributed  to  the  fact  that  the  profiles  seen  by  the 
main  landing  gear  were  rougher  In  the  latter  case.  Figure  7 shows  the 
Power  Spectral  Density  (PSD)  levels  of  each  line  of  survey  for  the 
Washington  National  runway.  A PSD  is  a measure  of  the  relative  roughness 
of  a runway  versus  frequency.  It  can  be  see*,  that  the  PSD  level  Is 
different  for  each  line  of  survey  which  accounts  for  the  change  In  the 
aircraft's  dynamic  response.  Figures  8 and  9 show  the  727-100  traversing 
a 1-cos  dip  headon  and  at  a 45®  angle  respectively.  In  this  case  the 
speed  W3S  50  fps  which  "tunes”  the  natural  pitching  frequency  (1  cps)  of 
the  727-100  to  this  1-cos  dip.  Hitting  the  1-cos  dip  at  an  angle  caused 
an  Increase  In  the  peak  P.S.  and  C.6.  acceleration  levels. 

In  was  necessary  to  try  to  simulate  different  aircraft  with  the 
computer  program  In  an  effort  to  check  the  program's  versatility. 
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TABLE  3 

SUMMARY  OF  SIMULA! lOMS 
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Boetng  727-100  Constant  Speed  Taxi  Simulation  over  the 
Washington  National  Profile  With  a Roll  Degree  of  Freedom 
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Figures  10  and  II  show  the  plotted  results  cr  the  C-9A  and  AMST  respectively 
hitting  the  1-cos  dip  at  a 45°  angle  at  a constant  speed  of  50  fps.  While 
the  C-9A  had  a relatively  high  response,  the  AMST,  which  Is  designed  to 
operate  off  of  rough  fields,  “absorbed"  the  dip  to  a large  degree.  This 
Indicates  that  the  computer  program  Is  calculating  relative  responses 
which  are  at  least  Intuitively  correct. 

Up  to  this  point  only  constant  speed  simulations  have  been 
discussed.  Figures  12,  13  and  14  show  the  plotted  results  of  the 
Boeing  727-100,  AMST  and  the  C-9A  respectively  taking  off  from  the 
Washington  National  runway  profile.  Takeoff  simulations  are  made  by 
starting  at  a near  zero  forward  velocity  and  accelerating  at  a constant 
thrust  until  rotation  speed  Is  reached,  then  the  simulation  is 
terminated.  Again  ft  can  be  seen  that  the  AMST  (idilch  is  designed 
for  rough  field  operation)  had  a very  low  dynamic  response,  even  though 
this  runway  Is  relatively  rough.  Figure  15  shows  the  plotted  PSD's 
of  Washington  National  Runway  36  and  of  two  typically  smooth  runways, 
one  at  Portland  Oregon  and  one  at  Dulles  International.  The  Washington 
National  Runway  Is  significantly  rougher. 

Experimental  data  was  available  for  comparison  with  the  727-100 
takeoff  slfiwlation.  Figure  16  shows  the  actual  time  history  plots  of 
vertical  acceleration  measured  on  a 727-100  taking  off  at  120,000  pounds 
gross  weight  from  Washington  National  Runway  36  on  December  11,  1972. 

Some  parameters  of  the  test  aircraft  were  unknown  such  as  strut  and 
tire  pressures  and  actual  Inertias.  So  exact  simulation  was  not 
possible.  However,  Table  4 shows  that  comparison  of  several  peak  values 
of  vertical  accelerations  at  the  P.S.  were  within  15  percent.  The 
comparisons  of  C.6.  vertical  accelerations  were  not  as  good.  In  the 
simulation  the  acceleration  levels  were  lower.  It  appears  that  main 
gear  strut  pressures  on  the  actual  airplane  were  lower  than  that 
simulated.  This  difference  would  cause  the  higher  response  in  the 
plunge  mode. 
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Figure  10.  C-9A  Traversing  a {1-cos)  dip  at  a 45®  angle 
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Figure  11.  AHST  Traversing  a (l-cos)  dip  at  a 45*  angle 
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Figure  12.  Boeing  727-100  Taking  Off  froei  Uashington  ffational  Airport 
With  £(^ie  Roll  Degree  of  Freedom  Included 
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Figure 


Washington  National  Airport  With  the 
Roll  Degree  of  Freedoju  Included 
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Figure  14.  C-9A  Taking  Off  from  Washington  National  Airport  With  the 
Roll  Degree  of  Freedom  Included 
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Figure  15.  PSD  of  Washington  National  Runway  36  and  two  Typically 
Smooth  Runways 
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Measured  Response  of  a Boeing  727-100  Takeoff  at 
Washington  National  Airport  Runway  36 
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TABLE  4 

COMPARISONS  OF  SIMULATED  AND  EXPERIMENTAL  DATA 


Experimental  Time 
(sec) 


P.S.  Vertical 
Acceleration 


C.G,  Vertical 
Acceleration 


Exp  Sim  Exp  Sim 

(g*s)  (g»s) 


8.0  0.9425  0.90  0.55 

16.3  1.305  1.12  0.80 

20.5  TAKEOFF 

Note;  All  measurements  arc  measured  from  peak  to  peak. 


0.35 

0.40 


The  remaining  simulations  made  were  with  the  inclusion  of  wing 
flexibility  in  the  simulations.  The  purpose  of  including  the  wing 
flexibility  was  to  see  if  there  was  a significant  change  in  P.S.  and 
C.G.  vertical  acceleration  response  when  the  wing  was  permitted  to 
bend  when  acted  on  by  a main  landing  gear  strut  force.  These  simulations 
were  made  on  the  C-9A  only  because  this  was  the  only  aircraft  for  which 
wing  flexibility  data  was  available.  Figure  17  shows  the  plotted  results 
of  C-9A  with  flexibile  wings  traversing  the  1-cos  dip  at  a 45®  angle. 

This  figure  can  be  compared  to  Figure  10  which  is  the  same  simulation 
without  flexible  wings.  By  superimposing  the  two  plots  it  was  determined 
after  T=17  seconds,  small  changes  in  vertical  acceleration  were 
appearing  in  both  the  P.S.  and  C.G.  responses.  Generally  the  higher 
accelerations  occurred  on  the  C-9A  simulation  with  flexible  wings. 

Also  there  was  a phase  lag.  By  the  end  of  the  run  the  rigid  wing  model 
lagged  the  flexible  wing  model  by  approximately  one  half  of  a cycle. 

Figure  18  shows  the  plotted  response  of  the  C-9A  with  flexible  wings 
during  a takeoff  simulation  from  Washington  National  .Runway  36.  Figure  14 
shows  the  same  simulation  without  flexible  wings.  Superposition  of  the 
two  plots  shows  little  change,  if  any,  in  the  airplane's  response. 
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Figure  17.  C-9A  with  Flexible  Wings  Taxiing  over  a {1-cos)  dip  at  a 
45"  angle 
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Figure  18.  C-9A  with  Flexible  Wings  Taking  Off  from  Washington 
National  Runway  36 
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SECTION  V 

SUMMARY  AND  CONCLUSIONS 

In  summary,  a mathematical  model  has  been  formulated  and  programmed 
for  a digital  computer  and  is  capable  of  simulating  most  flexible 
aircraft  traversing  an  unsymmetric  runway  profile  during  constant  speed 
taxi  or  takeoff-  Three  different  aircraft  have  been  simulated  and 
comparisons  have  been  made  with  experimental  data. 

Based  on  the  11  simulations  made,  the  following  conclusions  were 
drawn: 

1.  The  roll  degree  of  freedom  has  a significant  effect  on  the 
pilot's  station  and  center  of  gravity  vertical  acceleration  levels  if 
the  runway  profile  is  asymmetric.  The  degree  is  dependent  upon  how 
asymmetric  the  profile  is. 

2.  The  effect  of  wing  flexibility  on  T.S.  and  C.G.  vertical 
acceleration  r‘  .;^3nse  is  small  enough  to  be  neglected,  at  least  for 
the  airplane  simulated  (C-9A).  However,  with  the  addition  of  flexible 
wings,  it  now  becomes  an  easy  matter  to  expand  the  computer  program 

to  obtain  vertical  accelerations  (and  consequently  shears  and  moments) 
at  vital  wing  stations  such  as  the  wing  root  and  engine  and  stores 
pylons.  This  would  be  a natural  extension  of  the  study. 

3.  Comparison  of  the  simulated  aircraft  response  with  the 
limited  amount  of  available  test  data  was  satisfactory.  The  roughest 
parts  of  the  runway  were  identified  and,  as  in  the  test,  pilot  station 
acceleration  levels  exceeded  the  +0.4g  criterion.  If  exact  strut  and 
tire  pressures,  and  inertia's  were  known  for  the  test  aircraft,  the 
simulated  C.G.  response  may  have  more  closely  matched  the  experimental 
data. 


The  simulated  takeoff  took  an  additional  5 seconds  to  read 
rotation  speed.  It  is  assumed  that  the  actual  test  aircraft  we^  ht 
was  less  than  120,000  pounds,  because  several  runs  were  made  wi'  tout 
refueling  the  aircraft  after  each  run.  Therefore,  some  of  the  .uel 
had  been  burned  off.  The  fact  that  the  airplane  was  lighter  than 
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that  simulated  would  also  contribute  to  the  difference  In  C.G.  response. 
Also,  using  a 15“  flap  setting  changed  the  value  of  Cj^  and  resulted 
In  a shorter  takeoff  distance. 

4.  This  computer  program,  ''TAX2'',  appears  to  be  a very  efficient 
technique  for  lo rating  the  rough  areas  of  an  asymmetric  runway.  Using 
a CDC  6600  digital  computer,  a C-9A  takeoff  simulation  required 
70  seconds  of  central  processor  (CP)  computer  time,  which  Is  just 
30  seconds  over  real  time  for  this  simulation.  These  numbers  are 
typical  for  most  simulations. 


One  of  the  advantages  of  a program  of  this  type  Is  that  runway 
repairs  can  be  simulated  before  the  actual  repair  Is  made  In  order  to 
determine  the  minimum  amount  of  repair  required.  In  addition,  the 
effect  of  til?  proposed  repair  on  other  aircraft  can  be  determined 
before  the  repair  Is  made. 
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APPENDIX  A 

DEVELOPMENT  OF  EQUATIONS  OF  MOTION 

Development  of  equations  of  motion  using  Lagrange  equations.  All 
symbols  refer  to  Figure  A-1. 


'►  ◄ 
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dt  3q.  3q,  3q.  3^^ 


The  Kinetic  Energy  ?s: 


K.E.  - k 


The  Potential  Energy  is: 


P.E.  » + 


* ^ ♦ «Zc9  ■ ■•  ^C9 

* '“'‘sHL^^cg  * k ■ ^phl)‘ 

* '“''aiR^c,  - ^phr)^ 

* '^''sN^^cg  ■ ®®  ■ ^)  ♦ ■ Zrr)^ 


The  Dissipative  Energy  is: 


' * '“ilL'^cg  * 

* '■‘^KR'^cg  * Aa  - * i)‘ 


Now  Find  ^ 3KE 
^ 3^^ 


^ ^ % i i-  5ffi  . 


dt  3Z 


dt  91^, 


1.  Mi.  s >i,Zup  . jL  Mi  , I 0 

dt  32„p  ^ ’ 5?  ^ 


^ ^ ; i-  ^ . >xxi 

dt  aZj^  dt  3(j> 
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3K£ 


= 0 


Now  Find 


a(P.E) 

Jii 

3(P.E)  _ . 


'^'SHL^^cg 


AO  - Z„L 


C^) 


a(P-E) 

•:> 

'y.R 

3(P-E) 

'"ml 

9(P.E) 

3(P.E) 

3G 


^ W^cg"AO.Z„^  + C«) 

" V^cg 

= + Si  ■ '^SMR^^cg 

* ' ^PMR^ 

= + “m  ■ ^ML^^cg  " “ SiL  ' 

* ^M^SiL  ' 

* + Si  • '^SN^^cg  - BO  - Z^)  + - Zp^) 


= + ‘^SML^^^cg  AO  - Zn,l_  - C())) 
■*■  '^MRA^^cg  + Ao  - Z,^,p  + C({.) 


- KgjjB(Zcg  ’ Bo  - Z|^) 

iiMi  » . '^SML^^^cg  ■*■  AO  - - C(J)) 

^ ^MR^^^cg  + AO  - Z^P  + C*) 

Now  Find  3(D.E) 



. - 'HR<^eg  * - ^HR  + ci) 
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» - SiL^^cg  *0  ■ ^ ■ C.7) 
’"ml 

= . ^»cg  - M - ^|) 

ilMI  = , S„(Zcg  * "O  - ^R  * “) 

’^cg  • • . 

* SL^^cg  * ‘ 

- 'i.<^g  - • V 

= •*  W'^cg  ♦ *®  - ^R  * 

30 

* V<'cg  * k - «) 

- '^N^^cg  - 8®  - V 

ilp),  = ♦ C„RC(Zjg  » AS  - 2„„  ♦ W) 
^<^  * • • • 

■ V^'^cg  * *®  - Z„R  ■ "> 


Combine  Terms 


■ *^SML^^cg  ♦ A0  - - C<{i] 

■ ''sHR^^cg  * “ - ^R  * 8*3 

~ Sn^^9  * 80  - Zjj]  - H + L 

• ^ifMR)- 

■ SlR^^cg  ^ ^ " ^R 

_• 

' SlL^^cg  ■'■ 

* • 

" ^N^^cg  " • Zjj] 


Rewriting  we  have 


(1) 
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■ ^r(%)3  + 

- M + L 


where  the  terms  in  the  brackets  are  the  left,  right.and  nose  landing 
gear  strut  forces  respectively. 


Similarly 


' 'Vm^^MR  " ^PMR^  ' Hi 

W " S^V  “ - 'PN> 

■ W * 


xy 


^ = W<*«L>  * 'kf<’^L> 


■ "sHR*'**™*  ' W**nR* 


The  Forward  Translation  Equation  of  Motion  is  uncoupled  and 
expressed  as  follows: 


MX  =*  T - D,  - 


where: 


T * Thrust 


D,  = Aerodynamic  Drag 

O 


= Tire  Drag  (Total) 
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The  modal  method  will  be  used  to  express  the  aircraft's  flexibility 
as  follows: 


”i^'i  " ■ 2^<^iMiqi  - “i^  Mi^i 


where  i = the  ith  mode 

F.  = the  jth  force  input  into  the  system  (such  as  strut  force) 

J 

= the  itn  generalized  mass 

^l-j  = the  generalized  coordinate 

= the  modal  deflection  of  the  ith  mode  at  fuselage  station  j 
for  symmetric  modes  or  wing  station  j for  asymmetric  modes. 

u)^  = the  ith  mode  natural  frequency 

Z = Structural  damping  factor 

By  using  this  technique  the  displacements  X'^^  and  their  time 

derivatives  reflect  the  motion  of  the  bending  fuselage  and  wings  by 
N 

adding  the  2 q.?^^  (modal  displacement:-)  at  the  jth  (required 
i=l  ^ 

location). 

For  example; 

N P 

Total  Displacement  = Xj^jj^  + E + E 

1 —1  lc-*l 

. . ” . P . 

Total  Velocity  X'^^^  = X^,^  f E q^-^R  + E qj,S,^R 

i I K*”  I 


where; 


Term  1 = Displacements  of  the  rigid  body 

Term  2 = Displacements  due  to  the  symmetric  modes 

Term  3 = Displacements  due  to  the  asymmetric  modes 
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APPENDIX  C 

LISTING  OF  AIRPLANE  DATA 

- Boeing  727-100 

- McDonnell  Douglas  C-9A 

- AMST 
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APPENDIX  D 

LISTING  OF  RUNWAY  PROFILES 


1.  Figure  D-1  - Washington  National  Runway  36-Left,  Center,  and 
Right  hand  profiles* 

2.  Profile  data  listing  of  the  1-cos  dip  used  in  these  simulations! 
All  three  lines  of  profile  were  identical  except  that  they  occur  at 
different  times  corresponding  to  the  gear  locations. 
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